arXiv:1502.02048v2 [physics.chem-ph] 20 Feb 2015 


Molecular dissociation in presence of catalysts: Interpreting bond breaking as a 

quantum dynamical phase transition 

A. Ruderman 1 ’ 2 , A. D. Dente 3 , E. Santos 1,2 , and H. M. Pastawski 1 
1 Instituto de Fisica Enrique Gaviola (CONICET)-UNC and Facultad de Matematica, 

Astronomia y Fisica, Universidad Nacional de Cordoba, 5000, Cordoba, Argentina 
2 Institute of Theoretical Chemistry, Ulm University, D-89069 Ulm, Germany and 
3 INVAP S.E., 8403, San Carlos de Bariloche, Argentina 
(Dated: February 23, 2015) 

In this work we show that the molecular chemical bond formation and dissociation in presence of 
the d-band of a metal catalyst can be described as a Quantum Dynamical Phase Transition (QDPT). 

This agree with DFT calculations that predict sudden jumps in some observables as the molecule 
breaks. According to our model this phenomenon emerges because the catalyst provides for a non- 
Hermitian Hamiltonian. We show that when the molecule approaches the surface, as occurs in the 
Heyrovsky reaction of H 2 , the bonding H 2 orbital has a smooth crossover into a bonding molecular 
orbital built with the closest H orbital and the surface metal d-states. The same occurs for the 
antibonding state. Meanwhile, two resonances appear within the continuous spectrum of the ri¬ 
band which are associated with bonding and antibonding orbitals between the furthest H atom and 
the d-states at the second metallic layer. These move towards the band center where they collapse 
into a pure metallic resonance and an almost isolated H orbital. This phenomenon constitutes a 
striking example of the non-trivial physics enabled when one deals with non-Hermitian Hamiltonian 
beyond the usual wide band approximation. 

PACS numbers: 73.43.Nq, 34.35.+a 

INTRODUCTION 

When do two individual atoms become a molecule, or 
vice versa? This question is a fundamental one as much 
for Chemistry as for Physics. Certainly, it should involve 
a sort of discontinuity as in phase transition. In this 
context, P. W. Anderson hinted, in his well known arti¬ 
cle “More is different” [T], that a condition for a phase 
transition is the presence of infinite degrees of freedom. 

Sometimes, these are provided by the environment [2]. 

Many results in DFT calculations show an abrupt 
change in chemical bonds as the molecules approach to 
the surface of a catalyst and the interaction with the 
metal increases. In this context, a paradigmatic example 
of molecular formation and dissociation is the Heyrovsky 
reaction [3j, one of the steps of the Hydrogen evolution 
reaction at metallic electrodes: after the adsorption of 
a Hydrogen atom at the surface, a second proton ap¬ 
proaches and an electron is transferred from the metal. 

It is in this last step, when the Hydrogen molecule is 
formed and a discontinuity is hinted by DFT calculations 
as a jump in energy and spin polarization. This occurs 
at a critical distance of the farthest Hydrogen [3] . 

In Quantum Mechanics, a phase transition is recog¬ 
nized as a non-analytic behavior of an observable, typi¬ 
cally the ground state energy, as a function of a control 
parameter. This phenomenon, absent in a few atoms 
system, is an emergent property of the thermodynamic 
limit, i.e. when N, the number of atoms or degrees of 
freedom, tends to infinity. 0 ®. Such a limit is used to 
get the Fermi Golden Rule (FGR). There, one assigns 
an infinitesimal imaginary part —i 77 to each of the N 


involved energies and 77 is made zero only after taking 
N —i 00 [5]. Thus, an initial energy E 0 acquires a fi¬ 
nite energy uncertainty T associated with a decay rate 
2T/h. A straightforward way to account for this decay 
is to introduce an effective non-Hermitian Hamiltonian 
where some energies acquire an imaginary component, 
e.g. E 0 — > E 0 — ir. While this procedure dates back 
to E. Majorana (7], its deep physical implications only 
recently have become evident. This paper seeks to ratio¬ 
nalize the discontinuities found in DFT calculations at 
the light of the Quantum Dynamical Phase Transition 
(QDPT) concept [5]. 

A clear experimental evidence of a dynamical transi¬ 
tion showed up in NMR experiments in a 2(CH) 5 Fe crys¬ 
tal. There, the nuclear spins of the rare 13 C- 1 H dimers 
can perform Rabi oscillations. Beyond some critical crys¬ 
tal orientation, these spins are seen to abruptly decou¬ 
ple. This occurs when the interaction between the X H 
spin and the rest of the crystal becomes stronger than 
the 13 C- 1 H one 0 . This phenomenon is a QDPT and 
can also be interpreted as a particular case of the super- 
randiance phenomenon predicted by Dicke Hum. 

As a preliminary idea, let us explain how a tight- 
binding model for a homonuclear diatomic molecule A—B 
can show a non-analytical discontinuity. As in the spin 
case discussed above, this simple model would display 
the concepts relevant for the main discussion. Let SEq = 
2 Vab be the usual bonding-antibonding splitting. If 
atom B , has an interaction with a metallic band of width 
W SEq (wide band approximation) its energy acquires 
a finite lifetime H/2T due to the tunneling into the metal 
(r/2 <C W). This results in the effective non-Hermitian 


2 


molecular Hamiltonian [7], 


A MODEL FOR MOLECULE DISSOCIATION 


H c ff.= 


Ea —Vab 
—Vab E b — ir 


( 1 ) 


The eigenenergies are now complex numbers. For small 
r, the difference between their real parts is the splitting 
between the bonding and antibonding molecular levels, 

SE = \J[5Eq\ 2 — T 2 . This splitting is now weakened by 
the interaction with the substrate. Their corresponding 
imaginary parts T/2 are identical. However, SE has a 
non-analytic collapse when T reaches the critical value 
T c = SEq. Beyond this Exceptional Point (EP) [HH17] , 
the real parts become degenerate and the imaginary parts 
bifurcate. For big T values, the eigenvalue associated to 
B has an imaginary part T and thus remains strongly 
mixed with the substrate. The other one, associated to A, 
has an uncertainty proportional to |2Vab| 2 /F, indicating 
an almost isolated orbital [7 . A similar phenomenon un¬ 
derlays the Quantum Zeno Effect [18], i.e. when three or¬ 
bitals interact, a strong interaction between two of them 
weakens their interaction with the third. The detailed 
analytical and numerical solution of the above model was 
discussed in great detail in the context of QDPT by Dente 
et al. |19j . However, while very appealing for its simplic¬ 
ity, this picture can not he directly applied to a typical 
metallic catalyst because it is in a very different physi¬ 
cal regime. Indeed, the weak interaction with the wide 
sp-band does afford for a relevant role in the molecu¬ 
lar dissociation. Thus, our attention should turn to the 
strong interaction of the d-band with the molecule 20 ] . 
which has long been recognized as responsible for catal¬ 
ysis m- Yet the d-band has a width W d smaller than 
the molecular level splitting SEq, which prevents using 
the wide band approximation. Thus, we are back with 
the question of which is the magnitude that shows a non- 
analyticity that could be associated with a bond breaking 
in presence of a d-band? 

In this paper, we answer this question by showing that 
an actual analytical discontinuity appears if one includes 
a description of the metallic d-band with the right degree 
of detail. This requires to choose an appropriate com¬ 
bination between the molecular levels and the different 
metalic layers. Furthermore, the Anderson-Newns the¬ 
ory of adsorption mm is needed to describe molecular 
dissociation and electrocatalysis along with the inclusion 
of further works that extended this theory [23]. In our 
terms, while the molecule approaches to the surface, the 
farther A atom experiences a resonant through-bond cou¬ 
pling [M] with the 2nd layer of the metal. This interac¬ 
tion, mediated by the B atom and the first surface layer, 
manifests as two resonances inside the d-band. The tran¬ 
sition occurs when these resonances collapse at the center 
of the d-band, releasing the A atom. Meanwhile, the B 
atom hybridization with A is swapped into a H-metal 
bonding. 


We consider the Heyrovsky reaction, i.e. a H 2 molecule 
approaching perpendicularly to the metal surface. The 
molecule Hamiltonian is H$ = Ea |A) (A| + E b \B) (B\ — 
Vab (|A) (B I + | B) (A|). The degeneracy of the atomic 
energies Ea and E B is broken by the bonding inter¬ 
action Vab■ According to the standard wisdom [25J 
the d z 2 (top) or combination of the d xz and d yz (for 
hollow sites) are the only d-orbitals with a finite over¬ 
lap with the molecule. The orbital | B) interacts with 
the closest d-orbital combination, say 11 ), pointing along 
the connecting path through the binding energy Vo, 
Vsm = —Vo (| B) (1| + |1) (B |). The interaction energy 
Vo is roughly an exponential function of the molecule- 
substrate distance. Therefore, the complete Hamiltonian 
becomes H = H B + Hm + Vsm- 

Let us now focus on Hm which describes the d-band. 
Since the Newns pioneering work (21j . it is usually as¬ 
sumed that a semi-elliptical shape is a good approxima¬ 
tion for the Local Density of States (LDoS): 

N d (e) = x © [W d - 2e] x0[2e- W d ], 

(2) 

where 0 [a;] is the Heaviside function. Indeed, the va¬ 
lidity of this proposal for an actual metal, can be vi¬ 
sualized through the Lanczos transformation [26]. This 
is a unitary transformation which maps the actual 3 D 
substrate into an equivalent ID linear chain. Fig. [I] rep¬ 
resents this procedure for a 2D case. Starting from the 
| 1 ) surface d-orbital, a sequence of “collective orbitals” is 
build up mainly from the atomic orbitals at layers pro¬ 
gressively distant from the original surface orbital. Re¬ 
gardless of the precise original Hamiltonian Hm, in the 
Lanczos basis the metal Hamiltonian has identical 
“site” energies E n = E d and only nearest neighbor in¬ 
teractions V n ’s that account for the coupling among the 
Lanczos collective states [26]. Their rapid convergence to 
Voo = V = Wd /4 justifies the Newns proposal of a homo¬ 
geneous linear chain to describe the actual metal. The 
kets | n) are the collective d-orbitals, now the nth site of 
the metallic chain. Since in the surface of the metal the 
coupling between layers is weaker than in the bulk we 
assume |Vi| < |V 2 ] < \V\ and we set V n = V for n > 3. 

OO OO 

Hm = ^2E n \n) (n| + ^ -V n (|n) (n+ 1 | + |n + 1) (n|). 

n =1 n=l 

(3) 

With the purpose to define an optimal configuration 
for the molecular dissociation we will base our model in 
the Anderson-Newns theory and set up the Fermi energy 
level at E d , the center of the d-band, thus Ea = E B = 
E d - Thus the molecular bonding Vab yields a symmetric 
splitting around the center of the band. Additionally, we 
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2D Metal + Molecule Lanczos's Transformation 



FIG. 1. Effective non-Hermitian Hamiltonian: diatomic A-B 
molecule (in gray) in a configuration perpendicular to a 2D 
metal. The Lanczos unitary transformation combines orbitals 
at the same distance. The resulting tridiagonal Hamiltonian 
can be represented as a linear chain. A decimation procedure 
leads to a 4 x 4 non-Hermitian Hamiltonian. The same pro¬ 
cedure can be applied to the 3 D metal. Dot 1 represents a 
single atom for top interaction or a combination of them for 
hollow. 


will set the coupling elements V\jV = 0.8, V^/V = 0.9, 
and Vab/V = 2.5 which results fairly representative of 
various situations m- 

At this point it is necessary to list some of the main 
approximations implied by this model: 1) The fixed value 
of Vab neglects the variation of the distance between the 
atoms A and B. This does not affect our main results 
since the molecule breaking can be viewed as a compe¬ 
tence among interactions 2 ) As usual the atoms in the 
metal are considered fixed in the whole problem. How¬ 
ever, variations are minor in the Lanczos approach 3) 
We assume a null coupling between the metal and A the 
furthest atom. A residual exponentially small interac¬ 
tion would be completely masked by the through-bond 
interaction. 


MOLECULE DISSOCIATION AS A SPECTRAL 
BIFURCATION 

The solution of a linear chain model is better expressed 
in terms of the Retarded Green’s function matrix G = 
(el — H) “ 1 , whose divergences occur at the Hamiltonian 
eigenstates. For example, in absence of the metal we 
have: 




1 

e — Ea — LU(e) 


where 


Sa(e) = 


\Vab\ 2 
£ — Eb 


( 4 ) 

Clearly, Ea, the isolated A atom energy, is modified 
by the presence of B through the self-energy E^e), a 


real function accounting for the bonding and providing 
for the exact molecular eigenstates. As discussed before, 
one uses Ek = E^— ip (for k = A, B,n). This regulariza¬ 
tion energy 77 facilitates the study of the spectral density 
and the retarded charge dynamics, and whose physical 
origin can be traced back to small “environmental inter¬ 
actions”, as with the infinite sp-band states [ 55 ]. When 
the molecule is coupled with the metal in presence of the 
infinitesimal environment, one gets the retarded Green’s 
function [ 55 ] : 


Gaa(e) — 


£ - E a - 


| Vai 


£ — Eb — 


I Col 2 


e — Ei — 


\Vi\ 


M 

icr 


£ — £2 — ^tE(£) 


(5) 

Here, E(e) is the self-energy correction describing the 
bulk of the metal d-band in the Lanczos representation: 


S(e) = 


|U | 2 


£ - E d - 


From now on we set Ed = 
results [ 55 : : 


- = A(e) — ir(e) (6) 

E(e) 

0, thus, the solution of Eq. [ 6 ] 


s ( £ ) = £ —Y^--sgn{£)x ^ +isgn(y) x 

( 7 ) 

2 2 

£ — Tj £Tj / - 

with x = --- V 2 , y = -A and r = \J x 1 + y 2 . For 

the effective site at the second layer: 



G 22 (£) = 


£ + ip — E 2l{s) - 


M 

\V\‘ 


-E(£) 


with 


S 2 l(£) = 


|Ci 


£ T ip - 


|Co | 2 


£ + ip — 


\Vab\ 2 

£ + ip 


( 8 ) 


( 9 ) 


and similar equations for the other sites. In all these 
cases the LDoS is obtained from the imaginary part of 
the Green’s function TV, (e) = — y lim^o Im [Gu (e)]. 

The imaginary part of the self-energy accounts for the 
quantum diffusion of the electrons in the metallic sub¬ 
strate. Notice that the imaginary parts of the self-energy, 
T(e), and the continuum spectrum are already conse¬ 
quences of taking the thermodynamic limit of infinitely 
many sites in the chain. Otherwise E(e) would be a col¬ 
lection of divergences at discrete eigenenergies as E 2 l(e). 
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Within the d-band, the imaginary part survives the limit 
77 —>• 0 , indicating that each atomic orbitals merges into 
the metallic band [EE]- However, the mere existence of 
IV 0 does not warrant the QDPT. In this narrow band 
limit, a QDPT emerges as consequence of the specific 
non-linear dependence of T and A on £ that accounts for 
the different metal layers. 

Now, the important distinction respect to the intro¬ 
ductory example is that the self-energies acquire a non¬ 
linear dependence on e that contains all the wealth of the 
molecule-catalyst interaction. Finding the corresponding 
energy spectrum involves a 4th order polynomial on e 
with complex coefficients. A simple procedure is to find 
the eight complex roots of |1 /Gaa(e)| 2 - Half of them 
are non-physical as they are divergences for |Gaa(e)| 2 
but not for Gaa{z)- Thus, we evaluate the poles nu¬ 
merically. Once we obtain the solutions, we choose the 
physical ones, i.e. those whose imaginary component is 
negative, i.e. poles of the retarded Green’s function. 

In Fig. ?? we show the real part of the Green’s 
function poles. There, we observe two energies outside 
the d-band, which for Vo = 0 correspond to the bond¬ 
ing (|A) + \B))/y/2 (shorted as | AB)) and antibonding 
(|A) — \B))/y/2 (shorted as |(AU)*)) localized states of 
the lonely molecule (i.e. H 2 ) at tVab- When V 0 in¬ 
creases strongly, e.g. for V 0 = 3 V, these energies split fur¬ 
thermore and become a bonding and anti-bonding states 
between | B) and |1) orbitals, \B\) and |(H1)*), respec¬ 
tively. As in the quantum Zeno effect, increasing this 
interaction would dissociate the A atom from the rest of 
the system [T5] . This does not preclude a small amount 
of tunneling of A into the substrate’s second layer pass¬ 
ing through the orbitals B and 1. For intermediate Vo, 
this originates a through-bond interaction (211 that fa¬ 
vors the formation of a bonding state between |A) and 
12) , as well as an antibonding one(|A2) and |(A2)*) re¬ 
spectively). They are not localized states but resonances, 
as the electrons can be exchanged with the bulk. 

The above mentioned resonances, appear as poles of 
the Green’s function with a finite imaginary part ac¬ 
counting for their coupling with the metal, Fig. [3] How¬ 
ever, and this is the main point of this work, when Vo 
increases and reaches VJp, a quantum phase transition 
occurs and the state |A) becomes almost purely atomic. 
Simultaneously, state |2) recovers its purely metallic na¬ 
ture. At this transition the bonding and antibonding res¬ 
onances (identified by the real parts of the poles) collapse 
into a degenerate value. At this precise H 0 C the pole’s 
imaginary parts have a non-analytic bifurcation into a 
decreasing part, which accounts for a long lived atomic 
level, and an increasing uncertainty that represents the 
metallic delocalization of |2). This process can be inter¬ 
preted as a manifestation of the Quantum Zeno Effect 
meaning that when the interaction Vq between atom B 
and the metal increases, the interaction between A and B 
becomes weaker 0 HEj. In addition, the imaginary part 


Collapse of Resonances 



FIG. 2. Collapse of resonances in the d-band. The real part of 
the poles is plotted versus the molecule-substrate interaction, 
Vo- Small Vo values represent a far away molecule. For larger 
Vo’s the molecule is close to the metal. The gray area rep¬ 
resents the d-band, which goes from -2V to 2V. Outside the 
band the bonding molecular orbital | AB) smoothly becomes 
a bonding combination between the atom and the metal |£>1) 
(black curve). The same occurs for the antibonding state 
| (AB)*) that becomes | (£>1)*) (red curve). Poles within the d- 
band correspond to the bonding and antibonding resonances 
|A2) and |(A2)*} (green and blue curves respectively), after 
the transition they become an almost isolate | A) orbital and 
a 12) orbital strongly coupled to the metal. 

of the states out of the band remain zero for the whole 
range of Vo, due to the fact that they are localized states 
with an infinite lifetime. 

The detailed analysis of the spectral properties can be 
correlated with a study of the LDoS at different orbitals. 
Fig. a shows the results from such evaluation. In a) and 
b), we can see that for long distances (Vo — 0) the LDoS 
at | A) and | B) show a dominant presence outside the d- 
band. In contrast, c) and d) show that the sites 1 and 2 
of the metal do not have an appreciable participation at 
these energies. As the molecule approaches the surface 
and Vo increases, the LDoS at the orbital | A), shown in 
Fig. [4] a), looses its weight over the states outside the 
d-band. When Vo « V, we observe the emergence of 
two broad resonances inside the band accounting for the 
|A2) and |(A2)*) orbitals. Close to the non-analyticity 
point Vq , we observe that both resonances collapse into 
a single peak at e = Ea = 0. This is precisely the energy 
of an electron at the isolated orbital |A). 

An interesting complementary behavior is observed on 
the LDoS at |2). Fig. |4]d) shows that, as Vo increases, the 
two separate peaks, typical of a second layer J32], become 
close and almost collapse. They are still separated by an 
anti-resonance [24], i.e. a destructive interference with 
| A) which manifest as pole for Eq. [9j 

In Fig. [ 4 ] c) we observe that, as Vq increases, the 
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Bifurcation of T 



FIG. 3. Bifurcation of Decay Rates. The imaginary part of 
the poles is plotted versus the Vo molecule-substrate interac¬ 
tion. Small Vo values represent a far away molecule. Black 
and red lines over the abscissa account for the infinite lifetime 
of the localized states out of the band. At V<f the imaginary 
part of the resonances (blue and green curves) show a bifur¬ 
cation. One branch accounts for a the increasingly long lived 
atomic level, and the other branch describes the uncertainty 
of 12) that transforms it into the metallic delocalized band. 


a) Atom A b) Atom B 



e/v e/v 



e/v e/v 


FIG. 4. LDoS for different atoms and metallic sites, a) and 
b) A and B atomic orbitals, respectively, c) and d) 1 and 
2 effective metal orbitals, respectively. For atom A we ob¬ 
serve a decrease of the LDoS over the energies outside the 
band and an increment of its participation on the d-band res¬ 
onances. Atom B, instead, does not lose its participation over 
the states outside the d-band as Vo increases. The metallic 
orbital |1) loses its participation inside the band and gains 
presence over the localized states. The LDoS at orbital |2) 
shows almost no participation on the localized states outside 
the band. Consistently, the anti-resonance at the band cen¬ 
ters ensures no mixing with |A). 


first metal site starts loosing participation on the ri¬ 
band. Simultaneously, it increases its participation on 
the bonding and antibonding states localized outside the 
d-band. Accordingly, the | B) orbital maintains a sub¬ 
stantial weight in these localized states (Fig. @b) )• 
This accounts for the fact that after the phase transition, 
when atom A decouples from H, the out-of-band local¬ 
ized states become the bonding and antibonding |H1) and 
1(51)*). 


Resonances collapse inside the d- band 



FIG. 5. LDoS of A inside the d-band normalized to the 
local maximum. There we can see the two broad resonances 
collapse as Vo reaches Vo ■ This collapse describes the Quan¬ 
tum Dynamical Phase Transition. Furthermore, after Vf the 
LDoS becomes narrower as A becomes isolated. 

One might wonder how these prediction match with 
those of more realistic DFT calculation, such as the Ifo 
molecule interacting with a silver catalyst [2] . A first in¬ 
teresting effect, reproduced by DFT calculations and ob¬ 
served in the tight binding model, is the screening effect 
that the furthest atom suffers because of the presence ad¬ 
sorbed one. This result is observed in our LDoS on the 
bonding energies outside the d-band of Fig. [fj There, 
the participation of the | A) orbital is smaller than that 
of | B). However, the highly structured d-band masks the 
spectral branching that characterizes the QDPT evident 
in our tight-binding model, Fig. [5j after a renormaliza¬ 
tion to the LDoS maximum. Nevertheless, in a DFT, 
the hidden spectral bifurcation still constitutes the input 
for the full non-linear self-consistent calculation. Thus it 
triggers the discontinuity on some of the observables that 
show up. Indeed, while a Hydrogen atom approaches to 
another one adsorbed at the silver surface, a sudden jump 
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is observed in the total energy of the system. At the same 
time, the adsorbed atom also jumps to a position close 
to the approaching one, forming the H 2 molecule. Thus, 
the Heyrovsky reaction shows discontinuities whose roots 
may be assigned to the spectral discontinuities described 
by our QDPT model. 

CONCLUSIONS 

We analyzed the molecular dissociation in an Hetero¬ 
geneous Catalysis process under the framework of Quan¬ 
tum Dynamical Phase Transitions umm- As hinted 
by P. W. Anderson [3T|, the non-analyticity of the ob¬ 
servables is an emergent phenomenon enabled by an in¬ 
finite number of environmental degrees of freedom, in 
this case provided by the catalyst’s d-band. We first ob¬ 
serve a smooth crossover of the localized bonding and 
antibonding molecular states, which he outside the nar¬ 
row d-band, into a bonding and antibonding combination 
between the closest atomic orbital and the first layers 
of the metal d orbitals. By reducing the LCAO model 
to a non-Hermitian Hamiltonian where the imaginary 
parts have specific non-linear dependence on energy, we 
show that this system undergoes a collapse of resonances 
that provides the key to understand the dissociation phe¬ 
nomenon. More specifically, each of the resonances is 
formed from the bonding and antibonding interaction be¬ 
tween the furthest atom and a combination of d orbitals 
at the second layer of the metal. Before the molecule 
dissociation, both resonances are equivalently broadened 
by the rest of the metal. However, due to thw interac¬ 
tion with the surface they merge into a collective metallic 
molecular orbital centered in the second layer of the sub¬ 
strate and an isolated atomic orbital at the center of the 
d-band. 

In summary, we show that molecular dissociation con¬ 
stitute a striking example of the Quantum Dynamical 
Phase Transition, a simple but non-trivial phenomenon 
that only could emerge because we dealt with the ther¬ 
modynamic limit through non-Hermitian Hamiltonians. 
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